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emissions (during lifecycle and related to direct land use change), soil quality, water quantity and quality, 
and biodiversity. Suitable methods are selected and adapted based on an extensive literature review. The 


Key words: spatial variation in environmental impacts related to the spatial heterogeneity of the physical context is 
SEE S AN impacts assessed using Geographical Information System (GIS). The case study shows that there are large spa- 
GIS tial variations in environmental impacts of the introduction of bioenergy crops. Land use change (LUC) 
Miscanthus to sugar beet generally causes more negative environmental impacts than LUC to Miscanthus. LUC to 
Sugar beet Miscanthus could have positive environmental impacts in some areas. The most negative environmental 


impacts of a shift towards sugar beet and Miscanthus occur in the western wet pasture areas. The spatially 
combined results of the environmental impacts illustrate that there are several trade offs between envi- 
ronmental impacts: there are no areas were no negative environmental impacts occur. The assessment 
demonstrates a framework to identify areas with potential negative environmental impacts of bioen- 
ergy crop production and areas where bioenergy crop production have little negative or even positive 
environmental impacts. 
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1. Introduction 


In recent years, sustainability issues related to bioenergy have 
increasingly become a point of discussion in scientific, public and 
political arenas. At several levels, initiatives for sustainability cri- 
teria, codes of conduct and protocols have been developed to deal 
with this issue [1-5]. The steps towards sustainable biomass cer- 
tification system are well described by Lewandowski and Faaij [6] 
and van Dam et al. [7,8]. However, sustainable bioenergy criteria 
are currently hardly implemented in (inter-) national policies. In 
addition, initiatives regarding formulation of codes of conduct are 
mainly developed in a general and top down approach and do not 
account for the significance of regional variation in adoption and 
implementation. 

In the scientific arena, sustainability of bioenergy has been an 
important research objective. In Blottnitz and Curran [9] several 
studies regarding environmental impacts of bioethanol production 
chains have been reviewed. In most of the studies, energy balances 
and greenhouse gas (GHG) balances of production chains have been 
assessed. In others, such as Kaltschmitt et al. [10], Lewandowski 
and Heinz [11] Hamelink and van de Broek [12], Malca and Freire 
[13], Blottnitz and Curran [9] and Styles and Jones [14] more com- 
plete Life Cycle Assessments (LCA) have been performed. Van den 
Broek et al. [15], Smeets et al. [16,17] Rowe et al. [18] and Bor- 
jesson [19,20] are examples of studies on the sustainability of 
several bioenergy crops in specific countries that included a broad 
scope of environmental impacts and some even includes socio- 
economic impacts. The national scope of these latter assessments 
implies a relative aggregated level of analysis without differenti- 
ating for region specific conditions. Physical and socio-economic 
conditions often vary strongly between and within regions and 
determine the design of bioenergy chains: the chains need to be 
physically feasible in the geographical region, but also need to be 
compatible with the socio-economic context. In addition, the inter- 
actions between the bioenergy chain and the regional conditions 
have a strong influence on the actual environmental and socio- 
economic performance of bioenergy chains. A national approach 
could therefore give a distorted impression of the sustainability of 
bioenergy chains within regions. For that reason, the regional level 
is a more suitable level of analysis of the impact of biomass chains 
in order to safeguard sustainability of bioenergy chains for certi- 
fication purposes. This is illustrated by van Dam [21] who made 
a region specific assessment of impacts of bioenergy production 
in the Pampa region in Argentina. Although this study took dif- 
ferences in physical context into account, the impacts were not 
assessed spatially explicit. Regional averages on physical param- 
eters were used, e.g. for soil quality or water availability. As the 
spatial variation in determining parameters of impacts of biomass 
chains is significant, impacts would preferably assessed spatially 
explicitly. 

The aim of this paper is to develop a methodological framework 
to assess the environmental impacts of regional biomass chains 
in a spatially explicit and integrated way. The production of 
ethanol from Miscanthus (Miscanthus x Giganteus) and sugar beet 


(Beta vulgaris L.) in the North of the Netherlands (the province of 
Groningen, Friesland and Drenthe) is used as case study. The case 
area was selected because of the high data availability. This case 
study was also used to analyse the economic viability of regional 
biomass chains [22]. This study will explore the spatial variability 
of environmental impacts of bioenergy chains and builds on the 
previous analysis. 

The environmental impacts identified by several (inter-) 
national initiatives regarding sustainability criteria of bioenergy 
production were used as a starting point for a selection of rele- 
vant impacts and suitable indicators, as these criteria reflect the 
areas of concern for large scale bioenergy application [23-25]. The 
environmental impacts taken into account in this study are: GHG 
emissions (both during lifecycle and due to land use changes), soil 
quantity and quality, water use and water quality, and biodiver- 
sity. The selection and adaptation of methodologies to assess the 
environmental impacts is based on an extensive literature review 
on impact assessments. In this study, it is assumed that there is no 
competition with food production as the land for bioenergy crops 
is limited to the arable land that could become available by means 
of more efficient agricultural production. The developments in land 
availability for bioenergy crops in the Netherlands is based on the 
study of de Wit et al. [26]. 


2. Methods 


The level of spatial detail of the analysis depends on the avail- 
ability and resolution of the required input data. For several 
impacts, a postal code area was selected as the level of spa- 
tial analysis. For other impacts, a grid level of 0.1-1km? was 
selected. Because of the heterogeneity of the criteria, it is nei- 
ther possible nor appropriate to draw generic system boundaries 
and select a single functional unit. Therefore, system bound- 
aries and functional units were decided on for every individual 
indicator. 

In this study, it is assumed that only agricultural land (both 
arable land and pastures) is used for the cultivation of bioenergy 
crops. The effects of land use changes (LUC) for energy crop produc- 
tion were compared with the impacts of current agricultural land 
use. Nature areas were excluded beforehand. This is in line with the 
EU directive on the promotion of the use of energy from renewable 
sources [23]. The availability of agricultural land is therefore 
directly related to the intensification of agricultural land. As in 
the Netherlands the current state of agriculture is very advanced 
and consequently yield gaps are relatively small, only a small 
proportion of agricultural land is expected to become available for 
energy crops [22,26], Since current fertiliser and pesticide applica- 
tion levels in the Netherlands already meet legislative limits, it is 
assumed that the slight production growth of the remaining agri- 
cultural land are not a result of increased inputs but from improved 
management only. It is therefore assumed that this intensification 
will not result in significantly more environmental impact. Con- 
sequently, the environmental impacts related to the productivity 
increase of the remaining agricultural land were not included in 
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this assessment. The impacts of cultivation of Miscanthus and 
sugar beet were assessed for the entire agricultural area. This 
allows a straightforward interpretation of the analysis and pro- 
vides a rationale for the selection of most appropriate areas energy 
crop production. 

Only potential direct effects of implementation of bioenergy 
chains were assessed. Second order effects and feedback loops, e.g. 
the effect of changes in nutrient leaching on biodiversity, are hard 
to quantify and to allocate and are due to their complexity, beyond 
the scope of this study. Indirect effects and displacement effects 
due to LUC were tackled beforehand by assuming that only surplus 
agricultural land could be dedicated to bioenergy crops. Based on 
the first order results of the assessment of the spatial variability 
of environmental impacts combined with results of the study on 
the spatial variation in economic viability of bioenergy chains [22], 
areas can be identified where bioenergy crop production is eco- 
nomically viable and have the least negative or possibly positive 
environmental impacts. 

The environmental impact analysis of LUC to bioenergy crops 
on a spatially explicit and detailed level requires large amount of 
spatial data of soil characteristics, climate, current land use, yield, 
crop management, fertiliser and manure inputs and deposition lev- 
els, see, e.g. [16,21]. To deal with the large amounts of spatial data, 
a framework was developed to be able to calculate the impacts 
included in this study. 

The Miterra model used in our study is derived from the 
MITERRA-Europe model, which was developed to assess the effects 
and interactions of policies and measures in agriculture on nitrogen 
(N) losses on a regional level in EU-27 [27]. Miterra simulates the N 
and phosphorous (P) balance, emissions of NH3, N20, NOx, and CH4, 
leaching of N, NO3~ concentration in groundwater and changes in 
soil and biomass carbon stocks at a postal code level. For the appli- 
cation in this study, the model was adapted and input data and 
parameter values were made region specific. The main input data 
of the model are crop areas and livestock numbers, which were 
derived from the Geographical Information system Agrarian Busi- 
nesses database (GIAB) on postal code level. For the assessment of 
the impacts not included in the Miterra model, appropriate meth- 
ods were selected and adapted for spatial explicit application. In 
the following sections, the methods used to asses every individual 
impact are discussed. 


2.1. GHG emissions 


One of the main drivers of bioenergy is the reduction of GHG 
emissions. Therefore, the GHG emission related to bioenergy pro- 
duction is an important indicator to take into account. The GHG 
emission of the bioenergy chains was calculated by combining the 
GHG emissions related to direct LUC and the GHG emissions during 
all stages of the ethanol production chain. 


2.1.1. GHG during life cycle 

GHG is emitted during the lifecycle of biofuels and is related to 
cultivation of the energy crop (diesel for agricultural machinery, 
seed pesticide and fertiliser production and fertiliser application), 
storage and transport (heat for drying, diesel for transport, dry mat- 
ter loss) and processing (energy and chemical inputs for conversion 
process). The GHG emissions related to the production of bioen- 
ergy were calculated using a LCA approach which is often applied 
on bioenergy production chains. The ISO 14040-14049 guidelines 
articulate a preference for the use of system expansion in order to 
account for co products. As for legislative purposes energetic value 
based allocation is preferred in the EU directive on biofuels [28], 
this was applied in this study. The calculation of the GHG emissions 
was based on data provided by JEC [29] on the production of biofu- 
els. For cultivation parameters, universal data were substituted by 


region specific data provided in van der Hilst et al. [22]. In Annex 
1.1, the data used to calculate the primary energy requirements and 
the GHG balances are depicted. 


2.1.2. GHG emissions due to LUC 

GHG emissions due to LUC are caused by changes in soil car- 
bon stocks, above and below ground biomass and residues. When 
land is converted from one land use to another, carbon can accu- 
mulate (carbon sequestration) or diminish (carbon emissions). In 
addition, LUC causes changes in N20 emissions due to changes in 
fertiliser and manure application and drainage of organic soils. As it 
is assumed that there are no changes in livestock, livestock related 
N20 and methane emissions are assumed to remain constant when 
land is converted to energy crops. The Miterra model used here 
makes use of the methodology proposed in the IPCC guidelines to 
calculate the GHG emissions due to LUC [31]. In Annex 1.1, the 
parameters used to calculate changes in GHG emissions due to LUC 
are presented. 

Direct N2O emissions from managed soils related to different 
N sources (manure, grazing, mineral fertiliser, crop residues and 
cultivation of organic soils) and indirect N20 emissions due to N 
leaching and N deposition were calculated by Miterra using emis- 
sion factors of the IPCC guidelines [31,32]. The N surplus was 
calculated from the total N input (manure, mineral fertiliser, depo- 
sition and N fixation), the removal via crop harvesting and the N 
losses by gaseous emissions and runoff. Part of the N surplus is 
leached which results in direct N20 emissions [27]. In Annex 1.1 and 
1.3, the parameters to calculate the changes in N20 emissions are 
presented. 

Due to lower nutrient requirements, less manure will be applied 
on herbaceous energy crops. As no change in manure production 
is assumed, a manure surplus could be developed. This potential 
manure surplus is not considered here as it is outside the scope of 
this study. The annual changes in SOC, the change in carbon flux 
from organic soils and the changes in annual N20 emissions were 
combined to calculate the changes in annual GHG emission. 


2.2. Soil quality 


As soil quality is a broad and wide-ranging concept, it was nar- 
rowed down to indicators that are relevant for the region and for 
measuring the effect of bioenergy crop production on the soil. Soil 
organic matter (SOM) and soil erosion hazard are the indicators of 
soil quality that were selected in this study in line with the studies 
of Smeets et al. [33] and van Dam et al. [21,34]. 


2.2.1. Soil organic carbon 

As the SOC content is closely related to the amount of SOM, 
SOC is used as a proxy indicator. In addition to the benefits of car- 
bon sequestration in the soil (see Section 2.1.2), SOC is also linked 
to other important functions of the soil like water holding capac- 
ity, nutrient retention and soil structure [18,35]. SOC is therefore 
considered to be the most prominent indicator for soil quality [36]. 

SOC is affected by soil type, climate, past and present land use, 
soil and water management [37]. In the Netherlands, most soils 
have a high SOC content. Pastures usually contain more SOC than 
arable land due to continuous land cover, higher root biomass, high 
animal manure application rates and absence of tillage. Sandy soils 
have generally lower organic matter content than clay soils because 
of higher organic decomposition rates [37]. As wet conditions slow 
down decomposition rate, soils with high water tables contain gen- 
erally more SOC [37]. Peat lands store the largest amount of SOC but 
annual losses due to cultivation and drainage are high [37]. Changes 
in SOC were evaluated for the LUC related GHG. The methodology 
applied is discussed in Section 2.1.2. 
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Table 1 
parameters of the WEQ equation. 
Parameter Unit See Annex Ref. 
E Erosion Potential average annual soil loss tonha-!y 
I Soil erodibility index The potential annual soil loss per hectare for an open unsheltered tonha-!y Table A1.3 ap 
isolated levelled field with a smooth bare crustless surface. Related to 
texture, organic matter, calcium carbonate concentration 
K Soil surface roughness factor Reduction of the potential wind erosion due to ridges of the soil. factor Table A1.5 p 
Related to ridges and cloddiness made by field operations 
Cc Climate factor Local climatic erosivity. Related to wind speed, precipitation and factor Table A1.4 E 
temperature 
L Length of field Unsheltered distance across a field along the prevailing wind direction m 
V Vegetation factor Equivalent vegetative cover. Kind, amount and orientation of the tonha! y Table A1.5 d 


vegetation expressed in Small Grain equivalents (SGe) 


a USDA et al. (2002). 

b van Kerckhoven et al. (2009). 
€ Heijboer et al. (2002). 

d NRCS (2006). 


2.2.2. Soil erosion 

In areas prone to soil erosion, LUC may have a significant effect 
on the actual erosion risk. The main on-site problem caused by ero- 
sion on agricultural land is the loss of fertile top soil which leads 
to degradation of arable soils [37-40] and crop damage caused by 
abrasion or burial of seedlings or plants and the exposure and loss 
of seed [41]. In addition, the transport of minerals, organic matter, 
residues and pesticides could cause contamination of surround- 
ing surface water [39,42]. Furthermore, airborne particles due to 
wind erosion could affect human and animal health, machinery and 
infrastructure [39,40]. The main factors determining the actual ero- 
sion risk, are the soil characteristics (especially soil moisture and 
soil structure), vegetation (soil cover) and slope (in case of water 
erosion) [37,38]. It is therefore important to assess the effect of LUC 
or farm management on these factors in regions with high erodible 
soils and periods with high erosive rainfall or wind. 

Most studies regarding erosion focus on erosion caused by 
water, as this is the dominant form of erosion in most areas [38]. 
In the Netherlands, soil erosion due to runoff is mainly occurring in 
the sloping areas in the South-East [37,43]. The erosion hazard map 
of the Netherlands by Eppink and Spaan [43] shows no erosion risk 
due to rainfall runoff in the North of the Netherlands as this area 
has an outspoken flat morphology and the soils have a relatively 
high infiltration capacity. For this reason, rainfall related erosion 
is not included in this study. Wind erosion, however, does play a 
role in this region. Areas underlain by sandy and reclaimed cut-over 
peaty soils with frequent tillage on large fields without barriers are 
susceptible to wind erosion [37,40,43,44]. It is estimated that the 
vulnerable area of wind erosion in the in the North Netherlands 
amounts up to 60 kha. 

Several methods have been developed to model wind erosion for 
different temporal and spatial scales, functionalities (specific cir- 
cumstances) and impacts (soil loss, particle concentration). Ideally, 
soil erosion is continuously measured in the field or estimated by 
very exact wind erosion models fed with continues data on physical 
parameters. However, this is very time and capital intensive. In this 
study, the wind erosion equation (WEQ) method is applied. This is 
a relatively simple methods that requires less detailed input data 
and can be applied on a regional level as demonstrated by Van Ker- 
ckhoven [42]. The WEQ estimates the average soil loss due to wind 
erosion (in ton ha~! y~'), Eq. (1) and Table 1, give the parameters of 
the WEQ derived from USDA and NRCS [40], Morgen [45] and van 
Kerckhoven [42]. The complex relations between the parameters 
and the input data to calculate the soil loss are explained in Annex 
1:2. 


E= / (IKVCL) (1) 


When land is converted from conventional use to bioenergy 
crops, most factors included in the WEQ remain constant. The soil 
surface roughness factor could change slightly (420%) due to mod- 
ifications in tillage and planting practices. The main factor that will 
change is the equivalent vegetative cover. The erosion risk is deter- 
mined for every individual crop for every month of the year. To 
calculate the erosion risk for rotations, the risk of the rotation crops 
are combined in a weighted summation based on the portions of 
the crops in the rotation. For grass, it is assumed that pastures are 
reseeded every 10 years. For Miscanthus a lifetime of 20 years is 
assumed (1 year of establishment, 19 years of full crop cover, early 
spring harvesting). All parameters used for the calculations of wind 
erosion are depicted in Annex 1.2. In line with the WEQ guidelines, 
the most critical month has been selected based on calculations 
(see Fig. A2.3 and A2.4 Annex). 


2.3. Water quantity and water quality 


Water use and water quality are often addressed in proposed 
sustainability criteria for bioenergy [3,24,46]. As a relatively small 
proportion of the total water consumption for biofuel production is 
used during the processing stage (<10%)[47-50], only the water use 
related to crop production is taken into account. The assessment of 
the impact on water quality is limited to evaluation of the effect of 
change in use of agrochemicals like fertilisers on the water quality. 


2.3.1. Water use 

The change from current land use to energy crops may change 
the water balance of an area due to changes in evapotranspiration, 
runoff and percolation [33]. The amount of water lost through evap- 
otranspiration depends on crop type, growth stage, climate, soil 
characteristics growing period and agronomic practice [48,50-52]. 
As in the North of the Netherlands only a minor proportion of the 
agricultural land is irrigated (1% in wet years to 11% in dry years 
[53]), water withdrawal for irrigation is not incorporated in this 
study. 

Preferably, the impact of energy crops on fresh water availabil- 
ity for other functions is assessed on a water basin level [50]. This 
approach however, requires detailed knowledge and data about the 
hydrologic flows within a specific water basin. In this study, the 
changes in the water balance when current land use is converted 
to bioenergy crops was assessed on a spatial level of 1 km? grids. 
In addition, a comparison was made between the region specific 
water use efficiency (WUE) of the two bioenergy crops. 

To assess the potential water depletion due to the introduction 
of bioenergy crops, a simple water balance was made by compar- 
ing the evapotranspiration to the effective precipitation like done 
in the studies of Smeets and Faaij [17] and van Dam et al. [21] (see 
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equation in Appendix A1.3). The calculation method of water deple- 
tion is made time and spatially explicit by making use of spatial 
data on long term averages of monthly precipitation (1 km?) of the 
Royal Netherlands institute for Meteorology (KNMI). The effective 
precipitation (EP) is defined as the rainfall that is useful or usable 
in any phase of the crop production [54] and is derived from the 
actual rainfall making use of the USDA formula in the CROPWAT 
8.0 model [55]. The spatial distribution of monthly reference evap- 
otranspiration (ET)) was constructed on a 0.1 km? grid level for 
the North of the Netherlands by means of aggregation and spa- 
tial interpolation of daily ETg data of the KNMI weather stations in 
the region. Although the Penman-Monteith equation is the most 
widely applied method to calculate ETo, the KNMI uses the Makkin 
method for practical reasons [56,57]. 

The crop evapotranspiration coefficients (Kc factors) mainly 
depend on crop type, growth stage of crop and climate [58]. The 
crop specific decadal Kc factors for the Dutch situation were derived 
from the KNMI [59]. The spatial 0.1km? grid maps of monthly 
ETp were converted to land use specific monthly evapotranspira- 
tion maps, using the crop specific monthly Kc factors. The average 
monthly evapotranspiration of rotations were calculated using a 
weighted summation of the individual crops. 

As water deficits only occur during the summer in this area [59], 
the assessment of the change in total water deficit focussed on the 
months April-August. The parameters used for the calculation of 
the spatially explicit water deficit are depicted in Annex 1.3. 

The Water Use Efficiency (WUE) indicator is frequently applied 
in bioenergy related studies, like in Berndes et al. [48], Dorn- 
burg et al. [50], van Dam et al. [21], Fraiture and Berndes [47] 
and Smeets and Faaij [17]. In this study, the WUE is used as a 
second indicator of water consumption. It provides an indication 
about the water requirements per unit crop produced, whereas the 
water deficit methodology only provides figures for water use per 
hectare. In order to assess the spatial explicit WUE of bioenergy 
crops, knowledge on water availability (precipitation and ground 
water) and the effect of different growth limitation factors on evap- 
otranspiration rates should be considered. As here are too many 
uncertainties in these parameters, only aregional average WUE was 
calculated. 


2.3.2. Water quality 

The use of agricultural chemicals like pesticides and fertilisers 
may contaminate ground and surface water. A change in the use of 
these chemicals due to LUC in favour of energy crops is expected 
to have impact on the water quality. In this study, we focused on 
contamination of water due to fertiliser and manure use only, as 
these cause the most persistent problems in The Netherlands. An 
oversupply of N is likely to cause leaching of N to ground and sur- 
face water. The N surplus was calculated in the Miterra model as 
described in introduction section of the methods. Part of the sur- 
plus is denitrified and part will leach the ground and surface water. 
The amount of N leaching was calculated by multiplying the N sur- 
plus by a leaching fraction, which was based on the soil type, land 
use and groundwater level [60]. The NO3~ concentration in ground 
water was calculated by dividing the N leaching flux by the pre- 
cipitation surplus. For P there is no direct relation between the 
P accumulation in the soil and the concentration in ground and 
surface water. However, when P application levels exceed the P 
removal by crop uptake over longer periods of time, the risk of 
P leaching to water will increase. The P balance is calculated by 
Miterra making use of data on P inputs (manure, grazing, mineral 
fertiliser application and deposition) and P uptake by crops [32]. 
The input data required for the calculations of changes in NO3~ con- 
centration in groundwater and change in P balances are included 
in Annex 1.3. 


2.4. Biodiversity 


In most of the developed sustainability criteria for bioenergy, 
impacts on biodiversity have been identified as an important area 
of concern. LUC is a strong driver of changes in biodiversity [61-64]. 
Because of the loss, modification and fragmentation of habitats, the 
(indirect) expansion of agricultural land for energy crop production 
is perceived to be a major threat for biodiversity [65]. Biodiversity 
is also affected by LUC related depletion, degradation and pollu- 
tion of ecosystems and invasive species [64,66]. Changes in habitats 
due to energy crop production are most significant when natu- 
ral areas are converted to (intensive) agriculture areas [67]. Most 
bioenergy sustainability criteria deal with this issue by proposing 
process indicators, e.g. by referring to national regulations and by 
excluding protected areas and land identified as area with high 
biodiversity from bioenergy production [28]. In this study, these 
areas are excluded beforehand and only changes in agricultural 
land use are assessed. Currently, few guidelines are available about 
quantitative result indicators and methods to assess the impacts of 
energy crop production on biodiversity [46]. This is especially true 
for assessing agro-biodiversity on a regional level. 

Over the last decades, several indicator systems have been 
developed to assess changes in biodiversity. These indicator sys- 
tems vary to a great extent according to scale (global, national, 
regional or local), purpose (policy targets), and focal area of biodi- 
versity (species, genetic variation, population size or ecosystems). 
The impact of energy crop cultivation on biodiversity depends on 
both local scale effects (choice of crop, management intensity, veg- 
etation structure, substituted land use) and landscape scale effects 
(geographical location, scale and distribution of crops) [68]. The 
objective in this study is to quantify the impact of bioenergy crop 
production on agro-biodiversity on a local scale and identify the 
areas with a high risk on biodiversity loss. Two indicator systems, 
the High Nature Value (HNV) and the Mean Species Abundance 
(MSA) were selected for this assessment. 

HNV has become one of the main indicators to assess the impact 
on biodiversity of the 2007-2013 Rural Development Programs [69] 
and for the integration of environmental concerns into the Com- 
mon Agricultural Policy [67]. High Nature Value (HNV) farmland 
are certain types of farmland that, because of their characteristics, 
can be expected to support high levels of biodiversity of conser- 
vation concern [69]. JRC and EEA developed a method to identify 
HNV farmland [70]. Elbersen and van Eupen [71] adapted the HNV 
method for the Netherlands based on more detailed spatial data 
on HNV parameters. The HNV indicators used in this study are 
explained in Annex 1.4. Based on the characteristics of the HNV 
types and characteristics of the energy crops, the impact of a shift 
to bioenergy crop production on biodiversity in a HNV area was 
qualitatively assessed on a detailed geographical level. The risk 
level of biodiversity loss was identified based on the extent to 
which the new land use function could meet the habitat functions 
requirements of the important species identified by the HNV clus- 
ter indicators (see Table A1.9). The areas where the implementation 
of Miscanthus or sugar beet could cause a ‘high risk of biodiversity 
loss’, ‘considerable risk of a biodiversity loss’ and ‘both positive and 
negative effects on biodiversity’ were mapped. The HNV approach 
does not provide information about the impact of a shift towards 
bioenergy crops on biodiversity outside the HNV areas. In addition, 
the HNV approach does not measure the change in biodiversity 
quantitatively. 

The Mean Species Abundance (MSA) is a quantitative indicator 
for change in biodiversity. It does not reflect individual species 
responses but represents the average response of the total set 
of original species relative to their abundance in undisturbed 
ecosystems [72]. This indicator does not cover all aspects of the 
complex concept of biodiversity, but it can be used appropriately 
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Fig. 1. Overview of the bioethanol chain and the environmental impacts assessed, the indicators used and the models applied. The bottom boxes indicate the type of result 


that is generated from the assessments. 


to assess changes in biodiversity due to changes in land use for 
bioenergy crops. It was successfully applied in several global en 
regional studies concerning changes in biodiversity [50,73,74]. 
The MSA values for several land cover types are (included in 
Table A1.10 of the Annex) was linked to a map of current land 
use. However, as this land use map does not account for diversity 
within a land use type, management practices or locally available 
biodiversity, the MSA provides little differentiation for agricultural 
land in the Netherlands. For all areas outside the HNV areas, the 
MSA indicator is applied to assess the impact of LUC to bioenergy 
crops on biodiversity. 


2.5. Overview methods and integration of results 


Fig. 1 shows an overview of the bioenergy supply chain, the 
impacts that are assessed, the methods that are applied and the 
type of results that are produced in this study. 

For an integral picture, the GIS maps layers of the potential 
environmental impacts have been combined. As all impacts are 
measured in different units and scales, combining impacts requires 
standardisation of the GIS map layers. For all indicators, maximum 
standardisation was applied in which all negative impacts were 
translated to a scale from 0 to —10 and all positive effects to a value 
between 0 and +10. In this standardisation process, impacts equal to 


the impact of current land use were set to ‘0’, all increases in impacts 
compared to current land use were assumed to be negative and all 
decreases in impacts were assumed to be positive. In order to pre- 
vent a skewed distribution in standardised values due to extreme 
results of the indicators, the maxima were set on 2 times the stan- 
dard deviation from the average. In order to compare the results 
of Miscanthus and sugar beet, the same standardisation has been 
applied for the results of both crops. In this integrated assessment, 
it is assumed that all environmental impacts are equally important. 


3. Results 
3.1. GHG during lifecycle 


The GHG emissions over the production chain of bioethanol are 
substantial (see Fig. 2). The results found in this assessment are 
comparable with the figures found in previous studies [75-78]. 
Bioethanol production of Miscanthus causes fewer GHG emissions 
than first generation biofuels. As thermal energy requirements for 
2nd generation conversion are assumed to be fuelled by part of the 
feedstock, the process requires less fossil energy. Fig. 2 shows the 
relative importance of GHG emissions due to changes in SOC com- 
pared to the emission during the lifecycle. The uncertainty bars 
indicate the spatial variation in SOC change for a 20 year period. As 
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Fig. 2. GHG emissions due to SOC changes and during lifecycle of bioethanol production from Miscanthus and sugar beet. For comparison, the average GHG emission of 


petrol over the lifecycle is depicted. 


the emission reduction potential of a bioenergy production chain is 
highly related to the location of production this needs to be spatially 
assessed (see Fig. 3). 


3.2. GHG emissions due to LUC 


Fig. 4 shows the changes in GHG emissions due to LUC from 
current land use to Miscanthus and sugar beet. This includes the 
changes in carbon stocks (Fig. 3), the changes in carbon fluxes from 
organic soils and N20 emissions (Fig. A2.2 Annex). 

When grassland is converted to Miscanthus, a carbon loss due to 
land clearing (biomass carbon, up to 5.4 tonC ha~!) and soil prepa- 
ration prior to Miscanthus planting (soil carbon, 1 ton Cha~! on 
sandy soils, 2.5 ton Cha~! on clay soils) is compensated by relatively 
large above and below ground biomass stock (up to 17.5 ton Cha~!) 
that Miscanthus achieves compared to grass. When cropland is 
converted to Miscanthus, relatively large amounts of carbon are 
sequestered both in biomass and soil carbon (8.5tonCha~! on 
sandy soils, 21.5 ton ha~! on clay soils). For Miscanthus, the GHG 
emission reduction through carbon sequestration varies between 
0 to —78 kg GJethano!!. For sugar beet, the GHG emissions due to 
changes in carbon stocks vary between 0 and 54kgGJethanol!: 
These figures are in line with the studies of Foereid et al. [79] and 
Schneckenberger and Kuzakov [80] on sequestration potential of 
Miscanthus and with the study of Kuikman [35] on SOC in the 
Netherlands. 

LUC related N20 emissions contribute significantly to the GHG 
emissions. Conversion of current land use to Miscanthus gener- 
ally decreases N20 emissions. Especially when pastures (generally 
high fertiliser and manure application levels) are converted to 
Miscanthus, N20 emissions decrease significantly (41.7 ton CO2- 
eqha~! y—!). For Miscanthus, the GHG emission reduction though 
LUC related N20 emission reduction varies between 0 and —21 kg 
CO eq GJethano! ~. LUC to sugar beet cause an increase of N20 emis- 
sions in almost all areas, with highest increases on arable land. 
This is mainly caused by the higher fertiliser application levels and 


the increased emissions from crop residues (sugar beet leaves and 
crowns). For sugar beet ethanol, the LUC related N20 emissions con- 
tribute 0-15 kg GJetnano!” | to the GHG balance (see Fig. A2.2 Annex). 

Because LUC could cause an increase in one type of GHG but 
a decrease in another GHG, changes could balance each other out 
(for example: leaving the residues of sugar beet on the field will 
increase SOC but it will also increase N2O emissions from residues). 
In general, LUC towards Miscanthus will reduce GHG emissions, 
especially conversion of arable land. LUC to sugar beet will gen- 
erally cause a net increase in GHG emissions. The amount of GHG 
emissions is relatively high in areas now in use for pastures. For 
Miscanthus, the GHG emission reduction compared to current land 
use varies between 0 to —159kg COz-eq GJethanot | (—15 ton CO3- 
eqha~!). When arable land is converted to sugar beet, relatively 
small changes in GHG emissions occur. The GHG emissions due to 
LUC to sugar beet range from 0 to 148 kg GJethano!? (—1 ton CO2- 
eqha7!). 


3.3. Soil erosion 


The risk on erosion (in ton ha y~!) for current land use is shown 
in Fig. A2.5 in the Annex and the changes in the risk on erosion due 
to LUC are depicted in Fig. 5. The figures are based on the parameters 
for April which is the most sensitive period of the year for wind 
erosion due to low crop coverage, decreased soil moisture levels 
and relatively high wind velocities (See Annex 2.2). 

Although the overall risk on erosion is relatively small (max 
<6tonhay~!), this could still have severe negative effects. In the 
current situation, light sandy soils in the south east of the region 
(veenkolonién) are the areas most at risk for soil erosion. In the 
areas prone to soil erosion, relative large changes in erosion risk 
occur. When arable land is converted to Miscanthus, erosion risk is 
significantly reduced. Although it is assumed that grass is renewed 
every 10 years and Miscanthus every 20 years, risk on erosion 
increases when pastures are converted to Miscanthus. This is 
related to the planting date of grass which is generally seeded in 
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Fig. 3. (a) A SOC when current land use is converted to Miscanthus. (b) A SOC when current land use is converted to sugar beet. 


September. By the time the soil gets prone to erosion in spring, 
grass already covers the soil to a large extent. Miscanthus, however, 
needs to be planted in spring and has a slow crop development in 
the initial growth stage. Therefore, in the year of establishment the 
risk on erosion is relatively high. These findings are in line with 
Kort et al. [81]. The risk on erosion increases up to 3 ton ha~! y~! 
when sugar beet replaces crop rotations. When sugar beet replaces 
pastures, the risk on erosion can even increase up to 9tonha~! y~! 
on sandy soils. 


3.4. Water 
3.4.1. Water consumption 


In the current situation, evapotranspiration of crop rotation 
and pasture land exceed the effective precipitation in the months 


April to September (see Fig. A2.6 of the Annex). Therefore, a 
cumulative deficiency occurs during these months. As the effec- 
tive precipitation exceeds the evapotranspiration considerably 
during the rest of the year, all water shortages are eventually 
replenished. However, temporary shortages during the summer 
could cause damages to agricultural production or natural areas. 
The evapotranspiration of Miscanthus exceeds the evapotran- 
spiration of both pastures and rotations of arable crops from 
May to September. During July and August the evapotranspi- 
ration of sugar beet exceeds the average evapotranspiration of 
rotations and pastures. Therefore, LUC to bioenergy crops could 
cause higher water consumption rates during summer. As water 
tables are regulated artificially in the Netherlands, more water 
should be supplied to the areas with higher water require- 
ments. The developments in crop specific evapotranspiration 
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Fig. 4. (a) A GHG when current land use is converted to Miscanthus. (b) A GHG when current land use is converted to sugar beet. 


over time in relation to the effective precipitation is depicted in 
Fig. A2.6 of the Annex. 

The spatial variation in water deficits during the months 
April-August caused by evapotranspiration of current land use is 
depicted in Fig. A2.7 of the Annex. The largest water shortages occur 
in the western part, as precipitation is relatively low in this area 
compared to the rest of the region. The change in water shortages 
when current land use is converted to Miscanthus, are depicted in 
Fig. 6a. The spatial variation in change in water deficits due to the 
introduction of sugar beet is depicted in Fig. 6b. 

In addition to the spatial water balance, the regional WUE was 
measured. In Fig. A2.8 of the Annex, the WUE for Miscanthus and 
sugar beet is depicted. Most studies indicate that the WUE of C4 
crops is generally higher than the WUE of C3 crops [48,50]. However 


in this study, it was found that the WUE of Miscanthus (2.23 gkg~! 
water) is lower that for sugar beet (3.87gkg~! water). This can 
be explained by the fact that the competitive advantage of the C4 
photosynthesis pathway decreases significantly further away from 
the equator. In addition, as climatic conditions and agricultural 
management are optimal for sugar beet production, it achieves rel- 
atively high WUE. The WUE found for Miscanthus is in line with the 
2.1 gkg—! found by Clifton Brown [82] and the range of 1-9.5 gkg~! 
provided by Berndes [48]. The regional WUE of sugar beet is in line 
with the range of 2.55-3.83 gkg~! given by [50]. 


3.4.2. Water quality 
The Miterra model was used to calculate the nitrogen and phos- 
phorous balance in the soil. Based on the N balance and the water 
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Fig. 5. (a) A Wind erosion risk in kg soilha~! y~! when current land use is converted to Miscanthus. (b) A Wind erosion risk in kg soil ha~! y~! when current land use is 


converted to sugar beet. 


surplus, the NO3~ concentration in the ground water was calcu- 
lated. The maps of the N and P balance are included in Annex 
2.3. As N and P fertiliser application levels are closely linked, the 
spatial patterns of changes in N and P balances are quite similar. 
Because of the large differences between the current fertiliser and 
manure input on pasture and the fertiliser requirements of Mis- 
canthus, the N and P surplus will be significantly reduced when 
pasture area is converted to Miscanthus. Also in the areas where 
sugar beet substitutes grassland N and P surpluses are reduced. 
However, when sugar beet substitutes a rotation of crops, N and P 
surpluses increase. The high increases of N and P occur when rota- 
tions with high proportions of cereals are converted to sugar beet. 
This is mainly caused by the relatively high fertiliser application 
rates of sugar beet compared to other rotation crops. 


The differences in patterns between the N balance and NO37 
concentrations can be explained by the differences in leaching rates 
(lower for grassland, peat soils and lower precipitation surpluses). 
The NO3~ concentration in mg NO31~! is shown in Fig. 7. The cur- 
rent NO3~ levels are 72 mg!~! on average and 75% of the zip code 
districts exceed the standards of 50g1-! of the Nitrate directive 
[83]. In general, the NO3~ concentration decreases when current 
land use is converted to Miscanthus. The average change in NO37 
concentrations is —53 mg1~! and varies between 0 in areas with a 
large fraction of arable land and —187 mg1~! in zip code areas with 
a high proportion of grassland. When pasture area is converted to 
sugar beet, a decrease in NO3~ concentration (of max —75 mg1~') 
occurs. When arable land is converted to sugar beet, an average 
increase in NO3~ concentrations of 16.5 mg1~! is simulated. 


F. van der Hilst et al. / Renewable and Sustainable Energy Reviews 16 (2012) 2053-2069 2063 


p Pa 
ae 


Ses a 


Miscanthus 
A Water shortage (mm) 


E <--150 

EE -150--125 

GB -125 - -100 

E -100--75 B suit-up 
E -75--50 [Nature 
E] -50--25 J water 


Sugar Beet 
A water shortage (mm) 
E -< 
E o- 20 
E 20-40 
B -so E Buit-up 
= 60 - 80 = Nature 
HE 20-100 RGG 


0 10 20 40 Km 


Fig. 6. (a) A Water depletion during summer (mm) when current land use is converted to Miscanthus. (b) A Water depletion during summer (mm) when current land use 


is converted to sugar beet. 


3.5. Biodiversity 


In Fig. 8, the HNV areas were biodiversity is at risk when current 
land use is converted to bioenergy crops are mapped. It shows that 
in almost the entire pasture area in the north of the Netherlands 
which comprises high densities of important species (endangered 
bird species and large shares of European populations of meadow 
and/or wintering birds), biodiversity is at risk when current land 
use is converted to either Miscanthus or sugar beet. 

The HNV areas where biodiversity is at risk when converted to 
sugar beet (532 x 102 ha, of which 193 x 10? ha at ‘high risk’) are 
more extensive than for conversion to Miscanthus (473 10° ha, of 
which 91 x 10° ha at ‘high risk’). The areas where LUC to bioenergy 
crops is expected to cause a considerable risk of biodiversity 
loss, are considered to be areas where energy crops should be 


introduced with care. The areas with a high risk on biodiversity 
loss are considered to be no-go areas for bioenergy crops. In some 
areas (59 x 102 ha) the introduction of Miscanthus could have both 
a positive and a negative effect on biodiversity. These are mainly 
the heterogeneous areas and the extensive managed arable areas 
in the east of the region. 

For the areas not included in the HNV areas, the change in bio- 
diversity is expressed in the change in MSA value. This mapped in 
Fig. A1.11aand b of the Annex. The change in MSA value indicates an 
increase in biodiversity when pasture land is converted to Miscant- 
hus and especially when arable land is converted to Miscanthus. 
The MSA values indicate a decrease in biodiversity when pasture 
land is converted to sugar beet cultivation and a status quo in bio- 
diversity when sugar beet is cultivated in arable areas. The results 
show that sugar beet will induce a risk of biodiversity loss both 
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Fig. 7. (a) A NO3~ concentration when current land is converted to Miscanthus. (b) A NO3~ concentration when current land is converted to sugar beet. 


within and outside the HNV areas and that the only areas where no 
negative effects on biodiversity occur are the areas already in use 
for intensively managed arable crops. 


3.6. Integrated results 


In Fig. 9a and b the integrated results of all environmental 
impacts are spatially depicted for both Miscanthus and sugar beet. 
The results are expressed as an average score on a standardised 
scale from —10 to +10. In general, Miscanthus has better overall 
scores (mean =2.33) than sugar beet (mean=—0.85). The results 
discussed in previous sections show a wide distribution of scores 
for the individual impacts. As the integrated results show average 
scores close to 0, it illustrates that there are several trade-offs 


between impacts and that positive and negative scores balance 
each other out to some extent. 

There are no areas where only positive impacts occur when 
current land use is changed to sugar beet or Miscanthus. For both 
crops, the western wet pasture area and the peaty south east area 
of Drenthe appear to be the areas with most negative impacts. The 
conversion of pasture land to arable crops on (reclaimed) peat areas 
has negative impact on carbon stocks, risk on erosion and biodiver- 
sity. These negative impacts are relatively small for Miscanthus but 
quite severe for sugar beet. The impacts in these areas would even 
be more severe in case the pasture areas were not as intensively 
managed as they are today. Sugar beet has also a negative envi- 
ronmental performance in the eastern part of the region. This area 
is characterised by arable land on sandy soils. As sugar beet is a 
more intensively managed crop than average arable crop rotations, 
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Fig. 8. (a) Risk of biodiversity loss when current land use is converted to Miscanthus based on HNV cluster indicator. (b) Risk of biodiversity loss when current land use is 


converted to sugar beet based on HNV cluster indicator. 


negative impacts on the NO3~ concentration, P balance, soil erosion 
and also biodiversity are relatively large. For the same reasons, 
Miscanthus (extensively managed crop) scores relatively well in 
this area. Sugar beet has the least negative impacts in the northern 
pasture areas on clay soils and negative impacts are partly balanced 
out by reduced risk of nutrient leakage and water consumption. 
Comparing the integrated results of the environmental impacts 
assessed in this study with the results of the economic performance 
of sugar beet and Miscanthus in this region assessed in the study of 
Van der Hilst et al. [22], show that the areas with best environmen- 
tal performance are also the areas with the least cost of production 
(in €GJpiomass !)- These are the areas currently in use for intensive 


agriculture. However, in most of these areas bioenergy crops can- 
not compete with current land use as conventional crops have a 
better economic performance on these suitable soils. Based on the 
results of the environmental impacts and economic performance, 
an area of 24kha can be identified as most promising for Miscant- 
hus as in these areas it could compete with current land use, it has 
relatively low production costs and has the least negative and also 
positive environmental impacts. This well exceeds the proportion 
of agricultural land that could become available for bioenergy crop 
production in 2015 (9.6 kha arable land and 1.5 kha pasture land in 
the considered region) as indicated in the Refuel study [26]. As the 
study of Hilst et al. [22] showed that sugar beet cannot compete 
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Fig. 9. (a) Integrated result of environmental impact of a shift from current land use to Miscanthus based on standardised maps of individual environmental impacts. (b) 
Integrated result of environmental impact of a shift from current land use to sugar beet based on standardised maps of individual environmental impacts. 


with current land use and has relative high production costs in the 
whole agricultural area in the North of the Netherlands and this 
study shows mainly negative environmental impacts, it is assumed 
that sugar beet is not to be a promising bioenergy crop in this region. 


4. Discussion and conclusions 


Preferably, impacts would be assessed on most detailed spa- 
tial level. However, several input parameters of the Miterra model 
are only available on a postal code level. In addition, crop areas 
on a zip code level provides a proper indication of the crop mix 
in rotation, whereas a land cover recording of a single year on 
a detailed grid level ignores the effect of rotating crops. Further- 
more, as there are uncertainties in several key input parameters 


such as crop management, leakage fraction per soil type, water 
surplus etc., impact assessment on a detailed grid level could 
give a too optimistic impression on the accuracy of the results. 
There are some drawbacks of the Miterra model: The fertiliser- 
derived N20 emissions are calculated in accordance with the Tier 
1 approach of the IPCC. However, as indicated by Smeets et al. [84] 
and Leschen et al. [85] this neglects the variability of emissions 
due to differences in environmental conditions, management and 
crop type. The Miterra Model includes only regional averages of 
climate data and in the calculation of leaching of nutrients and 
decomposition of organic matter, water surplus of precipitation 
was included but water table levels were neglected. Despite these 
shortcomings, the Miterra model was able to combine several spa- 
tial datasets in a consistent way and therefore provides excellent 
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information on the relative differences in impacts on a spatially 
detailed level. 

As indicated by Hoefnagels et al. [76], Cherubini et al. [88] and 
van Dam et al. [89], GHG balances of biofuels are very sensitive 
for assumptions regarding system boundaries, allocation methods 
and input parameters. Therefore, the results of the GHG assessment 
should be interpreted with care. 

High soil losses due to wind erosion occur only incidentally 
when all relevant conditions are disadvantageous at the same time. 
By using an average monthly climate factor in the WEQ, no accu- 
rate indication is provided on how often erosion thresholds are 
met. Although the WEQ is not the most sophisticated approach to 
assess the risk on erosion, it does provide a method to discriminate 
between crop covers and to assess the spatial variability in the ero- 
sion risk on a regional level. It is recommended to use actual field 
measurements and more exact models to spatially assess erosion 
risk on a regional level in order to achieve more accurate figures on 
the exact amount of soil loss. 

The cumulative water deficit is assessed taking the 
spatial-temporal precipitation and evapotranspiration levels 
into account. It neglects however groundwater level, flow schemes 
of the water basin and the extensive water table regulation 
system. Therefore, the spatial water balance provides too little 
information to assess where actual droughts or water surpluses 
will occur. In order to assess the actual effect on water tables 
it is recommended that more advanced hydrologic models are 
applied. 

In this study the HNV and MSA indicators were applied to iden- 
tify areas with a high risk on biodiversity loss. The HNV approach 
enables to exclude the most vulnerable areas from a biodiversity 
perspective, which in line with the bioenergy sustainability cri- 
teria such as developed by RSB, RSPO, NTA8080. However, due 
to a lack of information on the response of species on the intro- 
duction of bioenergy crops, no quantitative assessment could be 
made. The MSA indicator on the other hand, enables to quantify 
the overall effect of a shift to bioenergy crops on biodiversity but 
because the MSA values does not account for diversity within a 
land use type, management practices or locally available biodiver- 
sity, the MSA provides little differentiation for agricultural land in 
the Netherlands. 

In order to assess the impact on biodiversity in a quantitative 
way, we conclude in line with Rowe [18] that more empirical data 
are required on the impact of specific bioenergy crops on the occur- 
rence of individual species. 

The calculations of the environmental impacts are all based on 
a very broad range of input parameters. As the results are affected 
by the assumptions made on climate, soil, agronomy and manage- 
ment applied, they should be interpreted with care. However, as 
the assumptions are consistently applied, it provides a good indi- 
cation the relative spatial differences in environmental impacts. As 
all impacts are related to ecosystem functionality, they are heav- 
ily interlinked with each other. In some cases, impacts reinforce 
each other and in other cases, trade offs between impacts occur. 
Although it is hard to prove causality and to quantify relations, 
it would be interesting to assess the spatial correlation between 
several impacts. 

All environmental impacts of Miscanthus and sugar beet cul- 
tivation have been integrated in a standardised map of overall 
environmental performance. This required an interpretation of the 
dose-response relationship: to what extend is the caused effect 
a negative or positive impact. Both the dose-response relation- 
ship and the thresholds for actual damage are expected to be 
spatially heterogeneous, time related, and, in some cases, con- 
tradictive. In addition, cause-effect relationships could be linear, 
quadratic, exponential, inversely proportional, etc. As there is 
too little knowledge available of the dose-response relationship, 


a simple approach of maximum linear standardisation with no 
threshold was applied in this study. More expertise is required 
on the translation of a change in value of an indicator to an 
actual impact. This requires not only knowledge on physical rela- 
tionships but also on the focal points of policy, the willingness 
to accept and the different points of view of the stakeholders 
involved. In the integration of environmental impacts all impacts 
were considered to be equally important, whereas ideally multiple 
perspectives from stakeholder on the hierarchy of environmen- 
tal impacts would be included in a spatially explicit Multi Criteria 
Analysis (MCA). 

This first explorative step could be further developed to 
an integrated assessment method by means of a stakeholder 
assessment and more advanced standardisation and weighting 
methods. 

The integrated assessment shows that there are no areas in 
the North of the Netherlands where only positive effects occur 
when bioenergy crops are introduced and that there are sev- 
eral trade-offs between the impacts. In general, sugar beet causes 
relatively many negative environmental impacts especially in pas- 
ture areas. In these areas, the LUC related GHG emissions are 
up to 148 kgGJethano!~! and the risk on soil erosion increases to 
9ton ha™!. Also, there is a high risk on biodiversity loss in these 
areas. On the other hand, a shift towards sugar beet could cause 
a decrease of 75 mg1-! in the NO3~ concentration of groundwa- 
ter and a decrease in the seasonal water deficits with 100mm. 
LUC to Miscanthus has less negative environmental impacts and 
in the areas currently in use as arable land many positive envi- 
ronmental impacts are expected to occur. In these areas, the LUC 
related GHG emissions are reduced by —159kg GJethano!?, the 
risk on soil erosion could be reduced with 4tonha~!, the NO3~ 
concentration is reduced with 53mglI~!, and the shift towards 
Miscanthus could have a positive contribution to biodiversity. On 
the other hand, at some locations the seasonal water depletion 
increases with 150mm arable land is converted to Miscanthus. 
For both crops, the western wet pasture areas appear to be 
the area with most negative impacts. It can be concluded that 
the spatial correlation between impacts and current land use is 
higher than between impacts and soil type. The combination of 
the spatial distribution of the environmental impacts and the 
economic performance of sugar beet and Miscanthus assessed 
in the study of Van der Hilst et al. [22] shows that the areas 
with best environmental performance are also the areas with the 
least cost of production (in €GJpiomass~!). Considering the large 
amount of negative impacts on areas currently in use as pastures, 
it could be recommended to exclude these areas for bioenergy 
crop production. In addition, as sugar beet causes relative many 
negative impacts, large scale introduction of this crop for bioen- 
ergy purposes should be carefully considered and implemented 
wisely. 

This assessment has resulted in understanding in which areas 
the introduction of bioenergy is likely to cause severe environmen- 
tal impacts and in which areas little negative and even positive 
effects are caused. This study well illustrates the high spatial vari- 
ability of environmental impacts of potential bioenergy chains and 
therefore it stresses the importance to assess the sustainability of 
bioenergy supply chains in an integrated and spatially explicit way. 
This study provides a set of methodologies to assess the spatial vari- 
ation of environmental impacts of bioenergy chains. Therefore, it 
contributes to identify promising locations for bioenergy produc- 
tion from an environmental point of view. Because of the generic 
characteristics, the methodologies and tools applied for this study 
can also be used for other case studies. It should be noted that 
this type of spatial explicit regional assessment on environmental 
impact of bioenergy chains cannot replace site specific environ- 
mental impact assessment. 
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